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Abstract 



(-H I We consider a degenerate partial differential equation arising in 

population dynamics, namely the porous medium equation with a 
bistable reaction term. We study its asymptotic behavior as a small 
parameter, related to the thickness of a diffuse interface, tends to zero. 
We prove the rapid formation of transition layers which then propa- 
gate. We prove the convergence to a sharp interface limit whose normal 

^ ' velocity, at each point, is that of the underlying degenerate travelling 

m ■ 
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1 Introduction 

In this paper we consider the rescaled porous medium equation with a 
bistable reaction term 

ut = eA(u™) + -f{u) innx (0, oo) 

(PI \^ = oua^x(0,oo) 

ou 

u{x, 0) = uo{x) in f] , 



^The first author is supported by the French Agence Nationale de la Recherche within 
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and study the sharp interface hmit as e ^^ 0. Here $7 is a smooth bounded 
domain in M^ {N > 2), z^ is the Euchdian unit normal vector exterior to dO, 
and m > 1. 

We assume that / is smooth, has exactly three zeros < a < 1 such that 

r(0)<0, /'(a)>0, /'(1)<0, (1.1) 

and that 

1 

mu"'~^f{u)du>0. (1.2) 



The above assumption implies that the speed of the underlying degenerate 
travelling wave is positive (see subsection 13. ip . so that the region enclosed 
by the limit interface is expanding (see below). This explains why the re- 
quirement (|1.2p is convenient for the study of invasion processes. 

As far as the initial data is concerned, we assume that < uq < M (for 
some M > a) is a (7^(17) function with compact support 

Suppuo := Cl{x € r2 : uo{x) > 0} CC ^ . 

Furthermore we define the initial interface Tq by 

Tq := {x ^ Q : uo{x) = a} , 

and suppose that Tq is a smooth hypersurface without boundary, such that, 
n being the Euclidian unit normal vector exterior to Eq, 

Eo CC and Vuo{x)j^O if x G Eq , (1.3) 

Uq > a in ^0 , uo < 0. in ^q , (1-4) 

where JIq denotes the region enclosed by Eq and Qq the region enclosed 
between dQ and Eg. 

Problem {P^) possesses a unique weak solution u^ as it is explained in 
Section [21 As e — )• 0, by formally neglecting the diffusion term, we see that, 
in the very early stage, the value of u'^ quickly becomes close to either 1 or 
in most part of 0, creating a steep interface (transition layers) between 
the regions {u^ ^ 1} and {u^ « 0}. Once such an interface develops, 
the diffusion term is large near the interface and comes to balance with the 
reaction term. As a result, the interface ceases rapid development and starts 
to propagate in a slower time scale. Therefore the limit solution u{x,t) will 
be a step function taking the value 1 on one side of the moving interface, 
and on the other side. 

We shall prove that this sharp interface limit, which we denote by Ej, 
obeys the law of motion 



iP' 



Vn = c* on Tt 



where Vn is the normal velocity of Tt in the exterior direction, and c* the pos- 
itive speed of the underlying travelling wave (see subsection 13. ip . Problem 
{P^) possesses a unique smooth solution on [0, T*""^) for some T"*"^' > 0. 
We denote this solution by F = Uo<t<T™''^ (r* x {t})- From now on, we fix 
< T < T'"*^^ and work on [0, T]. ~ 
We set 

QT:=nx (o,r), 

and, for each t € [0, T], we denote by 0.\ the region enclosed by the hyper- 
surface Ft, and by fij the region enclosed between dO, and Ft- We define 
a step function u{x, t) by 

(1 in of) 
u{x,t):={ fort€[0,T], (1.5) 

[ in of) 

which represents the formal asymptotic limit of n^ as e — > 0. 

Our main result, Theorem II. H describes both the emergence and the 
propagation of the layers. First, it gives the profile of the solution after a 
very short initial period: the solution u^ quickly becomes close to 1 or 0, 
except in a small neighborhood of the initial interface Fq, creating a steep 
transition layer around Fq {generation of interface). The time needed to 
develop such a transition layer, which we will denote by t^, is 0(e|lne|). 
Then the theorem states that the solution u^ remains close to the step 
function u on the time interval [t'^,T] {motion of interface). 

Theorem 1.1 (Generation and motion of interface) . Assume m > 2. Define 
/i as the derivative of f{u) at the unstable equilibrium u = a, that is 

fi = f'{a). 

Let r] € (O,min(o, 1 — a)) be arbitrary. Fix oq > arbitrarily small. 
Then there exist positive constants Eq and C such that, for all e £ (0,eo) 
and for all {x, t) such that t^ < t <T, where 



/i ^e| lne| , (1.6) 



we have 



u^(x,t) G < 



[l-2e,l + r]] if xEOf)\AAc,|i„,|(F 



[0, 1 + r/] if xen (1.7) 

[0,77] if xen['\Ma,{Ft), 



where J\fr{Ft) := {x € O : dist{x,Ft) < r} denotes the r-tubular neighbor- 
hood of Ft. 



Remark 1.2 (About the thickness of the interface). Since the construction 
of super-solutions is much more involved than that of sub-solutions, the 
statement (jl.7p is more accurate in fij than in Ql . More precisely, on 
the one hand, (|1.7|) shows that the convergence to 1 is uniform "inside the 
interface" except in 0{£\ lne|) tubular neighborhoods of the sharp interface 
limit; on the other hand, (|1.7p only shows that the convergence to is 
uniform "outside the interface" except in 0(1) tubular neighborhoods of 
the sharp interface limit. 

Remark 1.3 (About the assumption m>2). Note that the sub- and super- 
solutions we shall construct to study the motion of interface allow m > 1. 
Nevertheless, since we consider not well-prepared initial data, we need to 
quote a generation of interface result from [2] which is valid only for m >2 
(if 1 < m < 2 the partial differential equation is not only degenerate but 
also singular). When initial data have a "suitable shape", the restriction 
m>2 can be removed. 

As a direct consequence of Theorem 11.11 we have the following conver- 
gence result. 

Corollary 1.4 (Convergence). Assume m > 2. As e ^^ 0, the solution u'^ 
converges to u in Uo<t<T(f^f x {^})> where i = 0, 1. 

For the relevance of nonlinear diffusion in population dynamics models, 
we refer the reader to Gurney and Nisbet [H], Gurtin and Mac Gamy |1U| : 
density dependent diffusion is efficient to study the dynamics of a population 
which regulates its size below the carrying capacity set by the supply of 
nutrients. Since density dependent equations degenerate at points where 
ti = 0, a loss of regularity of solutions occurs and their support propagates 
at finite speed. 

Let us mention some earlier works on problems involving nonlinear dif- 
fusion that are related to ours. Feireisl [6] has studied the singular limit of 
(P^) in the whole space M-^, which allows to reduce the issue to the radially 
symmetric case. Hilhorst, Kersner, Logak and Mimura |llj have investi- 
gated the singular limit of the equation posed in a bounded domain of M , 
with f{u) of the Fisher-KPP type. Note that the authors in |11| assume 
the convexity of Oq which allows them to construct a single super-solution 
for both the generation and the motion of interface. Here, we dot not make 
such a geometric assumption. 

The organization of this work is as follows. In Section [21 we briefly recall 
known results concerning the well-posedness of Problem (/""). Section [3] is 
the body of the paper: we construct sub- and super-solutions to study the 
motion of the transition layers. Finally, we prove Theorem II. II in Section HI 



2 Comparison principle, well-posedness 

Since the diffusion term degenerates when u = a loss of regularity of 
solutions occurs. We define below a notion of weak solution for Problem 
(P^), which is very similar to the one proposed by Aronson, Crandall and 
Peletier [3] for the one dimensional problem with homogeneous Dirichlet 
boundary conditions. Concerning the initial data, we suppose here that 
uq £ L°°(r2) and uq > a.e. Note that in this subsection, and only in this 
subsection, we assume, without loss of generality, that e = 1, which yields 
the Problem 



(P) 



'uf = A('u'") + /(n) inf7x(0,oo) 
^ ^ - on 517 X (0, oo) 



du 
u{x,0) = uo{x) in f] . 

Definition 2.1. A function u : [0, oo) -^ L^(Q,) is a solution of Problem 
(P) if, for ah T > 0, 

(i) nGC([0,oo);Li(f]))nL°°(QT); 

(ii) for all tp G C'^{Qt) such that v? > and 7-— = on 50, it holds that 

ou 



n J Jqt J^ J Jqt 

(2.1) 

A sub-solution (respectively a super-solution) of Problem (P) is a function 
satisfying (i) and (ii) with equality replaced by < (respectively >). 

Theorem 2.2 (Existence and comparison principle). The following proper- 
ties hold. 

(i) Let u~ and u'^ he a sub-solution and a super- solution of Problem (P^) 
with initial data Uq and Uq respectively. 

If Uq <Uq a.e. then u~ < n^ in Qt ; 

(ii) Problem (P) has a unique solution u on [0, 00) and 

< n < max(l, ||uo||loo(q)) in Qt ] (2-2) 

(iii) ueC{Q^). 

The proof of the theorem above can be performed in the same lines as 
in [31 Theorem 5] (see also [13] and [5] for related results). The continuity 
of u^ follows from [S] . 

The following result turns out to be an essential tool when constructing 
smooth sub- and super-solutions of Problem (P^). 



Lemma 2.3. Let u be a continuous nonnegative function in Vt x [0, T]. 
Define Vl'l = {x e VL : u{x,t) > 0} and T^ = dn^ for all t G [0,r]. Suppose 
the family T := L)o<^t<T^t ^ {^} ^-^ sufficiently smooth and let v* he the 
outward normal vector on T^. Suppose moreover that 

(i) V(u"'') is continuous in 0, x [0, T] 

(ii) C[u] := ut - A(n") - f{u) = in {(x, t) G O x [0, T] : u{x, t) > 0} 



(iii) 



du: 



on dn*t, for all te [0,T] . 



Then u is a solution of Problem (P). Similarly a sub-solution (respectively 
a super-solution) of Problem (P) is a function satisfying (i) and (ii) — (iii) 
with equality replaced by < (respectively >). 

The proof of this result can be found in [11] . 



3 Motion of the transition layers 

3.1 Materials 

Underlying travelling waves. Hosono [12] has investigated travelUng 
wave solutions for the degenerate one dimensional equation 

Ut = (u'"),, + f{u) . 

He proved that there exists a unique travelling wave {c*,U), that the sign 
of the velocity c* is that of Jg u"^~^f{u) du, and that the profiles vary with 
the sign of the velocity. More precisely, for c* < 0, the front is smooth 
and U G C°°(M), whereas, for c* > 0, we only have (C/"*"!)' G L°°(K), 
but {U''^~^y ^ C(ffi). These different behaviors of the travelling waves are 
in contrast with the density independent diffusion models, where fronts are 
smooth whatever their velocities are (see [7]). 

In the present paper, the assumption (jl.2p implies that c* > 0. More 
precisely the following holds (see [12] for details). The travelling wave (c*, U) 
is the solution of the auxiliary problem 

' (C/™)" + c*[/' + /([/) =0 on(-c50,w) 

U{-oo) = 1 

U(0) = a 

(3.1) 
U' <0 on (-oo,w) 

(C/™)'(w) = 

U = on \oj, oo) , 



for some a; > 0. As z — >■ — oo, terms are exponentially decaying: 

max(l-[/(z),|f/'(z)|,|C/"(z)|) <Ce-^l^l for z < , (3.2) 

for some positive constants C and A. As z /^ to, we have 

lim (C/™-i)' (z) = -"^^—^e and lim(;7™-i)"(z) = _ (^ - 1) ^/(q) ^ 

(3.3) 
and U'{cj) € [— oo,0). Moreover, for a positive constant which we denote 
again by C, there holds 

\{U"'y'{z)\ < C\U'{z)\ = -CU'iz) for all z G (-oo,w) . (3.4) 

The cut-ofT signed distance function. Another classical ingredient in 
similar situations (see |14j or [8]) is a cut- ojf signed distance function d which 
we now define. Let d{-,t) be the signed distance function to Tt, namely 

dist(x,rt) for x G rjj , 



d(x,t) :=r "!"'!"' ^*? !""!%) (3.5) 



where dist(rE,ri) is the distance from x to the hypersurface Tt- We remark 
that d{-,t) = on Fi and that |V(i| = 1 in a neighborhood of the interface, 
say \Vd{x,t)\ = 1 if |d(x, t)| < 2do, for some do > 0. By reducing do if 
necessary we can assume that d is smooth in {(x, t) £ Q x [0, T] : \d{x, t)\ < 
3do} and that 

dist{Tt, on) > 3do for all t G [0, T] . (3.6) 

Next, let C{s) be a smooth increasing function on M such that 

' s if |s| < do 

C(s) = <^ -2do if s < -2do 
2do if s > 2do . 
We then define the cut-off signed distance function d by 

d{x,t) :=c(d{x,t)] . (3.7) 

Note that 

if \d{x,t)\<do then |Vd(x,t)| = l, (3.8) 

that d is constant (= 2do) in a neighborhood of dQ, and that the equation 
of motion {P^) yields 

if \d{x,t)\<do then dt{x,t) + c* = . (3.9) 

Moreover, there exists a constant C > such that 

\Vd{x,t)\ + \Ad{x,t)\<C for all (x,t) gQ^. (3.10) 



3.2 Construction of sub-solutions 

Equipped with the travelhng wave (c*, U) and the signed distance function 
d, we are looking for sub-solutions in the form 

n-(x,t):=(l-e)[/f^^^^^l±^i^^i^') =(l-e)C/(z,-(x,t)), (3.11) 



where 

z-ix,t):= '^^'^'^^f'''^P'\ (3.12) 

Lemma 3.1 (Sub-solutions). Let p > be arbitrary. Then, for e > small 
enough, u~ is a sub- solution for Problem (P^). 

Proof. In this proof (and only in this proof) we set u^ = u and z~ = z. 
Note that 

n*t ={xen: d{x,t) < -e\ln£\pe^ + euj}, (3.13) 

where Q* is defined as in Leinma [2.31 It follows that n = near the boundary 
dQ so that the Neumann boundary condition (Hi) in Lemma [2.31 is fulfilled. 
Since {U^y{uj) = we see that V(ti™) is continuous in Q. x [0, T]. Therefore, 
by virtue of Lemma 12.31 it is enough to prove that 

eC'[u] := eut - e^ /^{u"^) - f{u) < 

in {{x, t) : d{x, t) < —e\ Inejpe* -|- eu} = {(x, t) : z{x, t) < u}. 
By using straightforward computations we get 

eut = (1 - e) {dt + e\ Inejpe*) U'{z) 

e2A(n™) = (1 - e)™|V(i|2(C/'")"(z) + (1 - e)™eAd(C/'")'(z) , 

where z = z{x,t). Then using the ordinary differential equation ([/"^)" -|- 
c*U' + f{U) = 0, we see that 

eC'[u] = Ei + --- + Ei, 

with 

El := (l-e) [dt + c* - {1 - e)"'-^eAd{mU"'-'^){z) + e\ In e\pe^] U'{z) 

E2 := (1 - e)™(l - |Vd|2)(C/"^)"(z) 

Es:={{l-e)-{l-er){U"^nz) 

E^ := -f ((1 - e)Uiz)) + (1 - e)fiU{z)) . 

In the following we shall denote by C some positive constants which do not 
depend on e > small enough (and may change from place to place). 



We start with some observations on the term E4. Note that 

/ ((1 - e)u) - (1 - 6)/(n) = -euf'{9) + efiu) , (3.14) 

for some € ((1 — e)u,u). Hence 

I-E4I <Ce. 

Moreover, since /(I) = and /'(I) < 0, it follows from ()3.14p that, for u 
sufficiently close to 1, 

f{{l-e)u)-{l-e)f{u)>f]eu, (3.15) 

for some /3 > 0. Hence since [/(— 00) = 1, by choosing 7 ^ 1 we see that 

E4 < -l3eU{z) < --/3e for all z < -7 . (3.16) 

In the following we distinguish three cases, namely (|3.17p . (|3.19|) and (|3.2U|) . 

Assume that 

— e\ Inejpe* — ej < d{x,t) < —e\ lne|pe* + eu , (3-17) 

which in turn implies that —7 < z < uj. Since U' < on {—oo,uj) and 
U'{uj) € [—00, 0), it holds that U'{z) < —a, for some a > 0. If e > is small 
enough ([33]) shows that E2 = 0; from (jO]) we deduce that |£;3| < -CeU'{z); 
moreover we have |-E^4| < Ce. In view of ()3.9p . Ei reduces to 

Ei = {l- e) [-(1 - e)™-^eAd(mf/'"-i)(z) + e\ lne\pe'] U'{z) . (3.18) 

Since | - (1 - e)"*-ieA(i(mC/™-^)(z)| < Ce, inequality Ei < ye\lne\U'{z) 
holds. Collecting theses estimates we have 

C'[u\ <{ye\lne\-Ce)U'{z) + Ce 

< —jpae\ lne| + Ce 

<0, 

if e > is sufficiently small. 

Assume that 

- do < d{x,t) < -e\lne\pe^ - e^ , (3.19) 

which in turn implies that z < —7 so that (|3.16p implies £^4 < . Here 
again ([3^ shows that E2 = 0, from ([3^ we deduce that jE's] < -CeU'{z), 
and El reduces to (|3.18p . Hence we collect 

C^[u\ < U' (z) [-{I - e)eC + {I - e)p e\lne\ - Ce] <0, 



for e > small enough. 

Assume that 

-2do<d{x,t) < -do, (3.20) 

which in turn implies that, for e > small enough, z < —■^. In this range 
p.Sp and (j3.9p no longer apply but the exponential decay (|3.2p shows that 



\Ei\ + I-E2I + l^sl < Ce"^S. Last ^4 < -^/5e (see IKWh ) shows that, for 
e > small enough, C^[u] < 0. 

The lemma is proved. D 

3.3 Construction of super-solutions 

The construction of super-solutions is more involved: since we want them to 
be positive it is no longer possible to use the natural travelling wave (c*, U) 
which is compactly supported. Therefore we shall first consider slightly 
larger speeds c > c* which provide faster travelling wave solutions which 
tend to +00 in —00; then a small modification will provide us positive and 
more regular functions which are "nearly" travelling wave solutions. Before 
making this argument more precise, let us note that, as it will clearly appear 
below, the possibility of the above strategy follows from [12] . 



Let r] € (0,min(a, 1 — a)) be arbitrary. Let ao > be fixed. Let us 
recall that we have fixed < T < T"^"^ , where T™"^^ denotes the time when 
the first singularities occur in (P^)- Therefore we can select p > small 
enough so that the following holds. First the smooth solution (F^) of the 
free boundary problem 




on Tf 



exists at least on [0,T]. Secondly, if we denote by d^{x,t) the cut-off signed 
distance function associated with F'^ := Uo<t<T(rt x {t}) then, for all {x, t) £ 
Qt, 

d{x, t)>ao^ d^ix, t)>—. (3.21) 

Since c > c*, as explained in J12l Remark 3.1], there exists a faster 
travelling wave (c, V) which satisfies the same requirements as (c*, U) in the 
auxiliary problem (j3.ip . except that y(— 00) = -|-oo rather than C/(— 00) = 1. 
In particular, V is still compactly supported from one side. 

Next, for all n > 1, following the construction which comes before Propo- 
sition 4.1 in [12j (it consists in slightly modifying the above travelling wave 
(c, V)), we can consider (c, [/„) such that 

(i) Un satisfies the ordinary differential equation 

(C/„™)" + cUn' + f{Un) = on some (-00, Z„) , 

10 



where [/„' < holds 
(ii) Un is constant equal to some {6n] 
(iii) (t/„)-"i G CHR) 
together with Un{0) = a and C/„(— oo) : 



1 

m— 1 



> on [Zn, oo) 



-oo. Moreover (5„ — t- as n ^ oo, 



so that we can fix ng ^ 1 such that ((5„p)'"-i < r/. 

As a conclusion, if we denote Uno , ^no and Z„(j by W, 6 and Z we are 
now equipped with {c,W) such that W"^~^ G C^(]R) and 



(3.22) 



{W'^y + cW + f{W) = on(-oo,Z) 

W{—oo) = +00 
W^(0) = a 

W'<0 on(-oo,Z) 

1 
VF = (^"^-i < ry on [Z, oo) . 

We are now looking for super-solutions in the form 

dP{x,t) — e\ \ne\Ke^ 



uj{x,t) : 
In the sequel we set 



W 



z^{x,t) :-- 



dF{x,t) — e\ \iie\Ke* 



(3.23) 



(3.24) 



Remark 3.2 (The sub-domain S). We shall consider below a sub-domain S 
whose slice at time t, namely at := {x : {x,t) G S}, is the region enclosed 
between dVl and (more or less) T^. We shall prove that u^ is a super-solution 
in S. Thanks to (|4.4|) this will be sufficient for our purpose (see Section [l]). 

Denote by —9 the point where W{—9) = 1 + rj. For each < t < T, 
define the open set 

at := {x e Vt : (f{x,t) > e\\ne\Ke'' - e9} = {x : z^{x,t) > -9} , 

and the sub-domain 

S := Uo<t<r(fTt X {t}) . 

Note that the lateral boundary of T, is made of 9outS := dQ x (0,T) and 
5inS := Uo<i<T(st X {t}) where st denotes the smooth hypersurface 

st:= {x en-. (f{x,t) =e|lne|i^e* - e9} = {x : z^{x,t) = -9}. 

Lemma 3.3 (Super-solutions in S). Let r] G (0, min(a, 1 — a)) be arbitrary 
and let oq > be fixed. Then, for all K > 0, all e > small enough, uf is 
such that 



11 



(i) C^[ut] := {ut)t - £A(«)-) - -J{ut) > ^n S 

(ii) ^^„^^ ^ =0 on 5o„tS = afix(o,r) 

(iii) n+ = 1 + r/ on 5j„i; = Uo<t<T(st x {t}) . 

Proof. In this proof (and only in this proof) we put uf=u and z^ = z. Re- 
call that d^ is constant near the boundary dO, so that the Neumann boundary 
condition (ii) is fulfilled. Moreover the Dirichlet boundary condition (iii) is 
clear from the definition of st and the fact that W{—9) = 1 + 7/. 

Therefore it remains to prove that ££"^[1^] = eut — e^A(n"') — f{u) > 
in E = {{x,t) : z{x,t) > -6}. If z{x,t) > Z then C'iu] = £=[5^] > 0. 
We now assume that z{x,t) G {—6,Z), i.e. 

e\lne\Ke'' -e9 <d''{x,t) < e\lne\Ke^ + eZ . (3.25) 

Straightforward computations combined with the ordinary differential equa- 
tion [W"")" + cW + f{W) = yield 

eC'[u] = {dl + c)W' - e\ lne\Ke'W' - eAd'iW"^)' + (1 - \Vd^\^){W"')" . 

If e > is small enough, then ()3.25p combined with ()3.8p and ()3.9p — with d'^ 
and c playing the roles of d and c* — shows that the above equality reduces 
to 

eC'[u] = -eW' (I lne|Ke* + I\d^ {mW"'-^)) . 

Since W' < we have e£^[u] > 0, for e > small enough. 

The lemma is proved. D 



4 Proof of Theorem 11.11 

4.1 A generation of interface property 

We first state a result on the generation of interface. 

Lemma 4.1 (Generation of interface). Assume m > 2. Let ij > be 

arbitrary small. Then, for all x ^ Q, we have, for e > small enough, 

0<u%x,f) <l + i], (4.1) 

and there exists Mq > such that, for e > small enough, 

uoix) > a + Moe|lne| ^u''{x,f) >l-e (4.2) 

uo{x) < o-Moe|lne| ^u''{x,f) < e, (4.3) 

where f^ = //~^e|lne|. 

12 



Proof. We only give an outline since the arguments can be found in [T], [2]. 
Denote by Y{t; ^) the solution of the bistable ordinary differential equa- 
tion Yt = f{Y) on (0, cxd) supplemented with the initial condition Y{0;^) = 
^. Modulo a change of the time variable, we can use the sub- and super- 
solutions constructed in [2] to deduce that, for some C* > 0, for e > small 
enough. 



<u%x,f) < 



y(^;no(rE) + e2c*(e^*V^_l)^^ 



Next a straightforward modification of [U Lemma 3.9] — which specifies 
the instability of the equilibrium Y = a and the stability of the equilibria 
Y = 0, Y = 1 — shows that for e > small enough, for all ^ € (—1, 2), we 
have 

y(//-i|lne|;0 <! + ??, 

and that 

^>a + e\lne\ ^ Y{i^r^ \ In e\;0 > 1-e 
^ < a-e\lne\ ^ Y{i^r'^ \ In e\;^) <e. 

The combination of the above arguments completes the proof of the 
lemma. D 

4.2 Proof of the main theorem 

We are now in the position to prove Theorem 11.11 To that purpose we 
show that solutions are in between the propagation of interface sub- and 
super-solutions at time t^. 

Proof. Assume m > 2. Let rj G (0,min(a, 1 — a)) be arbitrary. First note 
that the comparison principle directly implies that inequality (|4.ip persists 
for later times, i.e. 

u%x, t) G [0, 1 -h 7]] for ah x € $7 and allt^ <t<T. (4.4) 

Since Vmq 7^ everywhere on Fq = {x € ^2 : uo{x) = a} and since Fq is 
a compact hypersurface, we can find a positive constant Mi such that 

if d{x,0) < — A/ie|lne| then ^0(2;) > a + Moe| lne| 
if (i(x,0) > Mie| lne| then ^0(3;) < a — Moe| lne| , 

with Mq the constant which appears in Lemma |4. 11 
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We first investigate the behavior "inside the interface" . In view of ()3.13p , 
we can choose p > large enough so that, for e > small enough, the sub- 
solution u~ defined in ()3.1ip is such that the set {x : u~{x,0) > 0} is 
included in {x : d{x,0) < —Mie\ lne|}. Therefore, from the correspondence 
(|4.5p . the estimate ()4.2p and the fact that u~{x, 0) < 1 — e, we deduce that, 
for all X G rj, 

uj{x,0) < u%x,t^). 

It follows from the comparison principle that 

u-{x,t-f) = (l-e)C/( -^^ ^ J <u''{x,t), (4.6) 

for aU (x,t) G X [t^,T]. We choose C » max{c* fi~^ , p e^ , A) so that 

-Cellnel +c*u~^e|lne| +e|lne|ne*"*^ C,, , 2,, , ,^ „, 

< --|lne| < --|lne|, (4.7) 



e - 2' ' - A' 

where A > is the constant appearing in (j3.2p . We take x G ^J \ 

■^CellnelC^t), i.e. 

d{x,t) <-Ce\\ne\, (4.8) 

and prove below that u^{x, t) > 1 — 2e, for t^ < t <T. Note that, for e > 
small enough, 

d{x,t-f) = d{x,t) + c*t'', (4.9) 

which, combined with ()4.6p and ()4.7p . implies 

u^(x,t) >(l-e)[/(-2|ine|) 

>(l-e)(l-C7e2) 

> 1 - 2e , 



where we have used ()3.2p . 

Last we investigate the behavior "outside the interface". Fix ao > 
arbitrarily small. For such ag > 0, we follow the strategy of subsection 13.31 
to construct super-solutions in S. More precisely define K = 2Mi, where 
Ml > is the constant that appears in ()4.5p : then construct uf{x,t) as in 
()3.23p . In view of Lemma 13.31 (ii). (in), the super-solution uf{x,t) and the 
solution u^{x,t + t^) satisfy a Neumann boundary condition on Oout^^ and — 
taking advantage of ()4.4p — are ordered on 9inS. Last we claim that (note 
that d'^{x, 0) = d{x, 0) since Fq = Fq) 

u^x, f) <W { -^^^ — = uf{x, 0) , (4.10) 
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for all X e ao = {x : d{x,0) > 2Mie|lne| - e9}. Indeed, ([13]) and (|0]l 

show that u^{x,t^) < e and the conclusion ()4.10p follows from the fact that 

1 
(Jm-i < ]/Y _ Hence the comparison principle yields 

u''{x,t)<W {^-^ ^ ^ J =u+(x,t-t^), (4.11) 

for all (x,t) e S with t" <t<T. We take x € 0J°^ \Mao{Tt), i.e. 

(i(x,t)>ao, (4.12) 

and prove below that u^{x,t) < r], for t^ < t <T. From (j3.2ip we deduce 
that (x,t) G S so that (|4.1ip applies. Note also that d'^{x,t — t^) = (F{x,t) + 
ct^. Therefore we infer from (|3.2ip that, for e > small enough, d^(x,t — 
t^) > ^ which in turn implies 

/^-2Mie'^e|lne|\ i 

uf{x,t-t^) <W ^ ^ =5^^ <r], 

since Vl^(+oo) = 6"^-^ . Conclusion follows from ()4.1ip . 

Theorem 1 1.1 1 is proved. D 
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